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Abstract. - It has been shown that a water drop can bounce persistently, when thrown on 
a super-hydrophobic substrate. We present here scaling arguments which allow us to predict 
the maximal deformation and the contact time of the drop. This approach is completed by a 
model describing the flow inside the drop, and by original experimental data. 



Introduction. - A liquid hall is a drop which remains quasi-spherical when brought into 
contact with a solid surface. Different examples of such objects have recently been described: 
let us quote pearl drops, which result from the extreme hydrophobicity of the solid ^L^SaJii^, 
liquid marbles, achieved by texturing the surface of the liquid [51ini[7], and (more classically) 
Leidenfrost drops, obtained by putting a small amount of volatile liquid on a very hot plate 
In such cases, quick transportation of tiny amounts of liquid becomes possible without 
any leak, which can be of great interest in microfluidics applications. At the same time 
these systems realize pure capillary "devices," and are worth being studied for their original 
properties. 

For example, a liquid ball (with a typical diameter 1 mm) impinging onto a solid substrate 
bounces off, as if it were a tennis ball hitting the ground pillUj. The rebounds are persistent 
(the restitution coefficient can be very high, of the order of 0.9), and observed in a large window 
of impact velocity. If the velocity is too small (typically smaller than a few centimeters per 
second), the drop gets stuck on the substrate. If too large (above around 1 m/s), the drop 
endures extreme deformation during the shock, and finally breaks into several pieces. In 
between, a more detailed analysis of the shock itself shows that the contact time of the 
drop with the substrate does not depend on the impact velocity V over a large interval of 
velocity [TT] . 

In this note, we discuss the maximal deformation of this sort of "water spring," and also 
the value of the contact time during the shock, in particular in the limit of elastic impact 
(small impact velocity). These predictions are compared with original experimental data. 
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Fig. 1 - Maximal deformation of a liquid ball impinging a solid at a small velocity (a: Ro ~ 0.4 mm, 
V = 0.094 m/s, W ~ 0.05); or at a larger velocity (b: Ro ~ 0.4 mm, V = 0.47 m/s, W ~ 1.). 



We shall mainly focus on the limit of small deformations where a drop of initial radius 
Ro bounces back with a high restitution coefficient of the order of 0.9. This corresponds to 
kinetic energies smaller than surface energies, i.e. to small Weber numbers, where the latter 
quantity is defined as 

W = pRoVy^. (1) 
This condition is achieved for small impact velocities, i.e. smaller than Vw'- 



For a millimetric water drop (7 ^ 72 mN/m, p = 10^ kg/m'^), Vw is of the order of 1 m/s. 
As seen in Fig. ^ the drop at its maximal deformation looks like a flattened sphere for small 
impact velocities {V < Vw), or a pancake for larger velocities {V ~ Vw). 

Scaling equations for a water spring (V <C Vw )■ ~ We start from the Euler's equation: 

Pjjj: = -^P + P9, (3) 

neglecting the effect of viscosity thanks to the high restitution coefficient. We consider the 
limit V <C Vw, in which the drop is deformed by a quantity X much smaller than the initial 
radius i?o (see Fig. ^a); the characteristic time of deformation scales as X/V. On the other 
hand, the Laplace pressure gradient scales as jX/Rq (for example, assuming an ellipsoidal 
form for the drop, the pressure jump at the equator or the apex scales as Ap ^ j/Ro{l±X/Ro) 
and changes over the length Rq). Thus, Eq. (O can be dimensionally written as 

pV^Rl ^ jX^ - pgRlX. (4) 

This equation also expresses the transfer of kinetic energy into the surface and gravitational 
terms associated with the drop deformation. If the velocity is zero, the drop is deformed by 
gravity by a quantity 5, as first discussed by Mahadevan and Pomeau 12; ■ which reads 

5^RIIk-\ (5) 



by use of the capillary length = yl/{P9), which is about 3 mm for water. 

We first consider the case where capillarity dominates gravity {X ^ 5). Then, we find 
from Eq. Q that the maximal deformation scales as the velocity: 

X^{pRlhf'v (6) 
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Fig. 2 - Analogy with a spring system: elongation x is plotted as a function of time t. The circle 
{t — 0) corresponds to the moment of impact and the square to the maximal deformation. The star 
thus indicates the moment of rebound, (a) general case, (b) small velocity limit, (c) large velocity 
limit. 

The contact time should be of the order of X/V, and thus can be written as 

r^{pRl/^f\ (7) 

Eq. Q, which yields a 1^-independent contact time of about one millisecond, indeed corre- 
sponds to the plateau observed experimentally 

Eq. lO allows us to make clear the limit of this regime. Since X must be larger than S, 
we find that the velocity must be larger than a characteristic velocity, 

V,r.{gRlK'f\ (8) 

which is about a few cm/s. If V approaches Vc, gravity cannot be neglected any more; we 
recast Eq. (0} as 

pV^Rl + 7(5^ ~ 7 (X - Sf (9) 

which expresses the energy conservation of an imaginary spring-mass system of initial velocity 
V and initial elongation S. 

Fig. 121 may help to have an intuitive picture for Eq. The elongation x of the imaginary 
spring is plotted as a function of time t, setting i = at the moment of impact. The left-hand 
side of Eq. is represented by a circle {x — 5, x — V a.t t — 0) while the right-hand side by 
a square {x — at t ^ t/'^)- The moment of taking off is marked by a star, and r is compared 

1/2 

in the plot with the period of the oscillator T ~ (^pR^/^^ 

Thus, we can graphically deduce the contact time t. In the small velocity limit V Vc 
{pV'^Rl <C jS'^), we see in Fig. Ejb that the contact time r is equal to the vibration period 
T. In contrast, in the interval commented above {Vc <^ V Vw) and with increase in V, 
T decreases from T and approaches the plateau value T/2. The complete variation of the 
contact time r as a function of the impact velocity V will be calculated more precisely in 
the next section. Note finally that, as stated above, the maximal deformation, in the limit 
of extremely small impact velocities, tends towards the constant 6, as seen from Eq. 0: we 
logically recover the static deformation. 

Local model. ~ The previous scaling arguments can be completed by considering local 
flows during the impact. We start from the incompressibility condition (V • v = 0), which 
reduces to 

1 drvr dvz „ 

+ 77^-0, 10 
r or oz 
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in the cylindrical coordinate system (r, 6, z) with vg = 0. Looking for a solution of the type 
Vr = Vr{r) and — Vz{z) satisfying appropriate boundary conditions, and introducing the 
equatorial radius R (see Fig. ^b), we find 

R R /i 1 \ 

= = -2— 2;, (11) 

setting 2 = at the substrate surface. The velocity potential (f) defined by v —Vcf) can thus 
be written as 

If the viscosity is neglected, the dynamics of </) is governed by Bernoulli's equation: 

p— + 2 + P + pgz = constant. (13) 

Here, the pressure p at the surface of the drop is given by the Laplace pressure, i.e. p — jC'ir, z) 
where the local curvature is denoted C. 

We evaluate Eq. (|13() at the apex (r = 0, z = ft.) and at the equator {r = R, z ~ h/2), in 
the limit of small deformation (a: ^ /Rq) for which h ~ 2Rq. (Here, x is defined as in Fig. 
^but not necessarily at its maximal deformation; the maximum magnitude of x is X.) To 
estimate the curvature difference AC = C {0, h) — C {R, h/2), we take the value in the static 
limit since we are in the regime of small velocities; the condition of constant pressure inside 
the drop can be integrated numerically, which gives AC ~ 6.8x/ Rq. In this way, we obtain 

-7pRlx/2 + 6.8-fx- pgRl = (14) 

Writing 6 ~ k^Rq/6.8 and ut'^ — 12.87/(7pi?o) (which gives as a plateau value for the contact 
time T/2 ~ 2.3-\/ pR^/j), we find a general solution: 

X — 5 = xq cos {wt + If) (15) 

The initial conditions lead to the relations cosip — 6/xo and sirn^s — V/ (loxo), from which 
we get xo — \/ S'^ + /ui^ and X = S + xq. The time of rebound r should be given by the 
equation, S = xq cos{ujT+ip), with tt < ujT+ip < 2tt. Hence, we find Lur+tp = 27r— arccos {S/xq) 
with if = a.TCCOs{S / Xo) {0 < f < tt) or 

T = 2tc I 1 arccos | — | | (16) 

where 

S/xo = (1 + {V/Vc)^) 

The analytical result Eq. is plotted in Fig. 13 and allows us to recover the predictions 
made via Fig. El the contact time decreases by a factor 2 when the impact velocity is increased 
from about Vc to Vw- 

Experiments. - We measured the contact time of bouncing drops at small impact ve- 
locities for millimetric water droplets hitting a super-hydrophobic substrate. Rebounds were 
sequentially recorded with a high-speed camera with a typical sampling time of lO^'* s for a 
contact time around 10~^ s, which gives a precision higher than 10 %. The results are plotted 
in Fig. ^ for two drop sizes. We indeed observed that the contact time significantly increases 
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Fig. 3 - Contact time r versus the impact velocity V in the elastic collision limit. The line is calculated 
from Eq. 11611 . The dotted line corresponds to the case of inverted gravity (see the text). 

Fig. 4 - Contact time of a water drops bouncing on a super-hydrophobic substrate as a function of 
impact velocity. Triangles correspond to a drop with _Ro ~ 0.6 mm, while circles to Ro ~ 0.4 mm. 
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Fig. 5 - Maximal deformation X of a water drop (_Ro ~ 0.4 mm) as a function of its impact velocity. 
Experimental data (circles) are compared with the scaling prediction (line) of Eq. (6) . 

Fig. 6 - Contact time of a water drop {Ro ^ 1mm) in the field of an inverted gravity. 



at low impact velocity (in the range below 10 cm/s). In addition, the ratio between the largest 
and the shortest times is found to be very close to 2, for both sizes. Finally, the velocity above 
which a plateau is observed increases with the drop size. 

All these observations are qualitatively in agreement with our predictions. A fully quantita- 
tive comparison would require a very accurate measurement of the drop radii (which determine 
the value of 14), and also taking into account the possibility for the drop to stick at very small 
velocity: for V ^ Vc, it is observed that the drop does not bounce off but remains stuck to 
the solid. In contrast, the comparison between the model and the data becomes quite precise 
at higher speeds: the plateau value is well described by Eq. {7)) with a numerical coefficient 
of 2.6±0.1 JJ, in good agreement with 2.3, the value obtained above. 

We also measured the maximal deformation X of the drop during the impact as a function 
of the impact velocity V (Fig. O. We found a linear behavior, in agreement with Eq. 

This first series of experiments confirms that the contact time of the bouncing drop deviates 
from its asymptotic value when we are in the regime of low impact velocity, i. e. in the linear 
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regime of deformation. We interpreted this effect as due to gravity, and tried to confirm 
this interpretation thanks to a second series of experiments. There, we studied the bouncing 
with gravity working in the opposite direction {inverted gravity). Drops were dropped from 
centimetric heights onto a super-hydrophobic substrate sUghtly tilted. After bouncing off, 
they hit at a height close to their maximum rebound height on a second plate of the same 
nature and inclined by the same angle. We recorded the contact time for rebounds on this 
second plate. The data are plotted in Fig. |H| as a function of the impact velocity. The 
error bars are there larger (because of the experimental resolution and the effects due to the 
vibrations induced by the first impact), but the data clearly indicate that the contact time 
increases with the impact velocity. This is in agreement with our theory: transforming g into 
—g (equivalently, S into —5) in Eq. (|16|l makes r increase with V as plotted in dotted line 
in Fig. 131 Interestingly, here, the contact time can vary between and r (between r and 2r 
in the previous experiments) which can be understood qualitatively using the construction 
suggested in Fig. [21 In the inverted gravity case, the whole sequence takes place below the 
line a; 0, so that we recover that the contact time should vary between and r. 

Perspectives. - Interesting perspectives can be given in the regime of higher impact 
velocities {V > Vw) where the drop is highly deformed and makes some kind of (transient) 
pancake of radius Rq + X X and thickness h (Fig. ^b). Since the drop crashes on the solid 
during a crashing time Rq/V, we expect the inertial term in the Euler's equation to be of the 
order of pV^ /Rq. The Laplace pressure scales as j/h at the equator and it changes over the 
length h. The pressure gradient thus scales as j/h'^. The gravity term can be neglected in 
this regime of large deformation, so that the Euler's equation can be cast dimensionally into 

pR',V'^jXyRl (17) 

via the conservation of the volume (i?Q ^ hR^). Thus, the maximal size of the pancake is 
deduced as 

X - RqW^I^ (18) 

where W is defined in Eq. QJ. X is indeed found to be larger than i?o for W > 1. We note 
that if the initial kinetic energy were transformed mainly into a surface term, namely, 

pRlV^ ^ ^X\ (19) 

we would find X - R^W^I'^ instead of Eq. (|TH|l. 

Eq. H17|l is in striking contrast with the case of a small deformation where the kinetic 
energy was found to be stored in surface (and gravitational) energy during the shock as in 
Eq. Namely, Eq. H17|l tells us that the initial energy is transferred primarily into other 
forms, identities of which (internal flow in the pancake, and possibly, the phonon energy of 
the substrate, etc.) remains to be clarified. 

In this regime, the contact time can be regarded as the dewetting time of the pancake. In 
this inertial case, the dewetting velocity Vd scales as yi/{pK) 0, which corresponds to the 
retraction speed of a hquid sheet ^Hl- Thus, the contact time is given by t X/Va^ which 
happens to result in Eq. 0. Indeed, the contact time was observed to be independent of 
V and to increase as i?Q^^ even for such high velocities IT!. This fact emphasizes that the 
drop is no longer in a linear-spring regime, for which Eq. H19|l would hold, and the contact 
time would be given by r ^ X/V . Note also that both X/Vd and X/V are certainly much 
longer than the crashing time Ra/V in the present case {W ^ 1), which justifies our previous 
estimation of the inertial term in Eq. H17|l . 
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Conclusion. - In the limit of small deformations (i.e. for small impact velocities), a 
liquid ball thrown on a solid behaves as a quasi-ideal spring. This can be understood as a 
conventional spring-mass system with a stiffness given by surface tension and a mass given 
by that of the ball; the deformation of the small ball during the impact linearly depends on 
the impact velocity and the contact time scales as the period of this spring-mass system, as 
observed with high speed photographs. The contact time was found to increase (typically by 
a factor of 2) at small impact velocity, which can be interpreted as the result of the weight of 
the ball. This was confirmed by achieving a similar experiment in an inverted field of gravity, 
which indeed leads to a decrease of the contact time at small velocities. 

In this system, the effect of viscosity could be ignored. This might be physically due to 
the absence of a contact line in a situation of non- wetting; most of viscous dissipation usually 
takes place near the contact. It would be interesting to see how our views would be modified 
when the viscosity of the liquid consisting the ball is increased. Extensions to similar (but 
different) systems such as gel balls 14^ or drops of surfactant solution would also be worth 
studying. 
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